Statistics of resonance states in a weakly open chaotic cavity 
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In this letter, we demonstrate that a non-Hermitian Random Matrix description can account for both spectral 
and spatial statistics of resonance states in a weakly open chaotic wave system with continuously distributed 
losses. More specifically, the statistics of resonance states in an open 2D chaotic microwave cavity are investi- 
gated by solving the Maxwell equations with lossy boundaries subject to Ohmic dissipation. We successfully 
compare the statistics of its complex-valued resonance states and associated widths with analytical predictions 
based on a non-Hermitian effective Hamiltonian model defined by a finite number of fictitious open channels. 

PACS numbers: 05.45.Mt,05.60.Gg,03.65.Ad 



Examples of waves in enclosures can be found in as di- 
verse contexts as room acoustics, guided optics, vibrations of 
structures, etc. While these situations imply different physical 
mechanisms, they all belong to the domain of Wave Chaos, 
which reveals universal features of cavities with non trivial 
geometries ||XJ] . As long as these systems can be considered 
as closed, their generic spectral and spatial properties are cur- 
rently well described through the theory of Hermitian random 
matrices from a statistical point of view. Nevertheless, no re- 
alistic system is truly closed, thus calling for a description of 
the coupling mechanisms to the environment. Physically, the 
latter include bulk absorption, leads or waveguides, as well 
as dissipative or radiative boundaries. In the domain of Wave 
Chaos, open systems are actively investigated both from ex- 
perimental and theoretical points of view (see Refs J2, 0] for 
reviews). Among the domains concerned by experimental 
studies, one can cite: microwave cavities fl, optical micro- 
cavities Jfl, and elastodynamics JUS]. 

To analyze these open chaotic systems, the scattering ap- 
proach was found to be a powerful theory J^. In this frame- 
work, a system composed of N resonances is analyzed in 
terms of a N x N non-Hermitian random matrix, the so-called 
effective Hamiltonian: 



(1) 



where its Hermitian part H corresponds to the Hamiltonian of 
the closed system and its anti-Hermitian part — h VV"f models 
the coupling to the environment. More precisely, the open- 
ness is introduced by means of the N x M coupling matrix 
V whose elements connect the n = 1, . . . , N states to the 
j = 1,...,M scattering channels jj]]. The non-Hermiticity 
of HeB yields a set of complex eigenvalues {E n } associated to 
two distinct sets of eigenvectors called left { (ip n | } and right 
{\4>n)} eigenvectors: 



FIG. 1 . Real (on the left) and imaginary (on the right) component of 
the 500th resonance state obtained by means of the Finite Element 
Method. 



form a bi-orthogonal and complete set which can be normal- 
ized by: (iJ n \lpn) = Snm- 

While field statistics of open chaotic systems have been 
systematically studied for a given energy and considering 
the energy as a continuous parameter (see J3l and references 
therein), statistics of resonance states i.e. left and right eigen- 
vectors of H c ff, for TRS systems are less understood. For 
these systems, the impact of the openness is to turn real eigen- 
functions into complex internal wavefunctions associated to 
resonances. This complexness being uniquely related to the 
presence of currents inside the system 0,01 (after a phase ro- 
tation leading to independent real and imaginary components 
J3l). In order to quantify the presence of currents, one can 
use the complexness parameter introduced by Lobkis and 
Weaver ifioll as the ratio of the variance of the imaginary and 
real parts of the nth resonance state. 

Making use of right eigenvectors: 



E,(Re[^]) 2 



(3) 



Kefi\i>n) = £n\lpn) > (VU?4if = (i>n\£n , 
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where the eigenvalue £ n = E„ — iT n /2 gives respectively the 
energy E n and the resonance width T n of the ?ith resonance. 
The left and right eigenvectors, which describe the resonance 
states respectively for systems with positive or negative gain, 



where corresponds to the ith component of the right eigen- 
vector (we note that q\ can also be defined using left eigen- 
vectors). This parameter has recently gained attention both 
experimentally M fTlTl and theoretically 1 12, 12] in particular 
due to its relationship to the resonance width. A linear rela- 
tionship between q n and r„ was first noticed by Barthelemy 
et al. Illll analyzing hundreds of resonance states of a 2D 
chaotic microwave cavity at room temperature. This result 
was then confirmed using the effective Hamiltonian formal- 
ism in the limit M 3> 1, relevant in the experiment Il2ll . 
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A linear relationship between q n and T n was also verified in 
an elastodynamics experiment for a given resonance when a 
spatially extended coupling is varied jg]. Lately the complex- 
ness parameter was investigated at arbitrary M by means of its 
probability distribution in the regime of weak coupling 1 1311 . 
There, it was shown that the average value of q\ is directly 
proportional to the variance of r„, which constitutes the nat- 
ural measure of the fluctuations of the widths. 

It is the aim of the present work to confront our theoretical 
predictions to numerical solutions of the Maxwell equations in 
a 2D chaotic microwave cavity with lossy boundaries subject 
to Ohmic dissipation. After a brief introduction to the theo- 
retical model 11311 . we will describe the cavity we numerically 
investigate. Then we will discuss the statistical results we ob- 
tained concerning the complex wavefunctions and the widths 
of its resonance states and will compare them to the theory. 

In the regime of weak coupling, which was shown 11411 to 
correspond to the condition ^/var(T) <C A, where A is the 
mean level spacing, the anti-Hermitian part of W c g is small 
compared to the Hermitian part and the perturbation theory 
can be applied. As the eigenvalues are nondegenerate with the 
probability 1 due to the linear level repulsion at small spacings 
111], one gets directly expressions of the spectral widths and the 
complexness parameters (|3): 
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np 
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where Y np = Ylj=i ^n'^p an d tne {^n} correspond to the 
energies of the closed system. We would like to stress that g 2 
also gives information about the non-orthogonality of eigen- 
functions, with important implications in various physical sit- 
uations. For example, the non-orthogonality induces an en- 
hancement of the line width of a lasing mode 11511 . influences 
branching ratios of nuclear cross sections 1 kJ and is also in- 
vestigated in open quantum maps 11171 Il8ll . The distribution 



of the energies {E n } corresponds to the eigenvalue distribu- 
tion of the Gaussian Orthogonal Ensemble with N — > oo |Q]]. 
The coupling amplitudes are chosen to be real Gaussian ran- 
dom variables with zero mean and co variance (VfiVfy — 
a 2 S nrn S^ k [19]. To obtain an analytical expression of the 
distribution, we use the fact that T np = Y^jLi ^n^p can 
be viewed as a scalar product between M-dimensional vec- 
tors and then can be expressed using polar variables fl9ll : 
Tnp = y/TnTp cos 9 np . The complexness parameter is now 
given by: 
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p^n 



(5) 



where A np = T p cos 2 6 np 



Considering that the widths are 
given by a sum of M squared independent Gaussian random 
variables the distribution of the rescaled widths 7 = T /a 2 
(with (7) = M) is given by a \ 2 distribution with M degrees 
of freedom: 



X 2 m(7) 
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Xm=i which is in - 



and a further calculation shows that A satisfies the Porter- 
Thomas distribution Q: P(A/a 2 ) 
dependent of M. The distribution of the rescaled parame- 
ter X n ee £ql is defined by Vf° E {X) = (5(X - X n )), 
where the statistical averages are performed over the energies 
{E n }, the resonances widths {T n } and the {A,,}. By making 
use of group integral methods 11511 it reads: 
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where strong mode-to-mode fluctuations clearly appear, em- 
bodied in the power law tail 1/X 2 of the distribution. 

In the following, this prediction will be compared with the 
numerical solutions of the Maxwell equations in a 2D chaotic 
microwave cavity. The chaotic cavity we consider has a quar- 
ter of a stadium shape with a radius of R = 1 m and a length 
of I = 2 m (see Fig. Q]). In order to reduce the bouncing 
ball modes between the two parallel sides, an oblique cut is 
performed on one side and a movable perfectly reflecting half 
disk of diameter d = 0.3 m is placed on the opposite side. The 
absorbing boundary condition is imposed on an adjustable 
part of the upper arc of circle of the cavity. For the TM po- 
larization, the electromagnetic field (E, H) inside the cavity is 
uniquely characterized by the single component of the electric 
field: ip(r) = E z (r), where r = (x, y). The field component 
is solution of the Helmholtz equation: 



(8) 



where cl> is a complex angular frequency and c the light ve- 
locity. Starting from the ideal closed system where the field 
satisfies the Dirichlet boundary condition i.e. an infinite con- 
ductivity along the whole contour of the cavity: ip \c = 0, 
the openness is introduced through a finite conductivity a c on 
a length l a b s of the contour. To first order, the electric field ip 
along the lossy contour is then given by: 
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(9) 



where a c is the effective conductivity of the contour, and n 
is the unit normal vector directed toward the interior of the 

states, the Finite Element 



conductor 111 , 

As we examine resonance 



Method 12111 reveals to be very efficient to solve the time in- 
dependent wave equation (0. Using the commercial software 
Comsol™, we are able to obtain the complex eigenvalues 
u> n = oj n — i(, n /2 and the resonance states ijj n of the first 800 
resonances. To increase the statistical samples at our disposal, 
ensemble averages were performed by sliding the half disk re- 
flector along the largest side of the cavity (ensuring a constant 
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area and perimeter of the cavity) with 7 different positions 
well enough separated to produce completely statistically in- 
dependent spectra leading to samples of 2100 resonances for 
each numerical distribution. To compare the numerical results 
from the electromagnetic cavity with the theoretical model, 
the usual correspondence between (E n , T„) and (uj^, io n ( n ) is 
performed. 

In the 2D chaotic cavity we study, the number of channels 
M is related to the effective absorbing length l aos along the 
boundary. This number can be evaluated by using the Sabine's 
law of reverberation known in room acoustics [12J. This law is 
formally equivalent to the so-called Weisskopf's estimate for 
the level width (well-known in Nuclear Physics 1U9IP provided 
the effective number of channels be related to the absorbing 
perimeter through the intuitive relationship : 



M 



labs 

A/2' 



(10) 



where A = 2-kc/uj is the wavelength. Clearly, as this estimate 
is frequency dependent, to compare numerical distributions 
to our theoretical predictions, one must use samples of res- 
onances for which the parameter M is approximately a con- 
stant. This can be achieved by considering resonances within 
high-frequency intervals. Indeed, as the cumulated number of 
levels grows like cu 2 and M like to (according to dTol)), the rel- 
ative variation of M within a sequence of A7V adjacent levels 
around the Nth level is given by AM/M = ±AN/N. In 
practice, we considered intervals of 100 adjacent resonances 
above the 300th resonance. To explore various numbers of 
channels, namely two frequency ranges were considered : 
from the 300th to the 400th resonances (AM/M ~ 14%), and 
from the 700th to the 800th resonance (AM/M ~ 7%), and 
three different lengths of the absorbing part of the perimeter : 
labs = 7r/18, labs =7r/6 and l aos — tt/2 (in unit of radius). It 
is also important to note that the theoretical distribution of the 
complexness parameter (0 being established in the regime of 
weak coupling, the conductivity of the absorbing part of the 
perimeter was chosen such that ^/var(r)/A < 10~ 2 . 

Fig. |2] shows the distributions of the spectral widths z = 
7/ (7} of the cavity compared to the \ 2 distribution Four 
different sets of values of the absorbing lengths, conductivities 
and frequency ranges are used (see caption). For each interval 
of frequency, the number M of channels used for the compar- 
ison corresponds to the nearest integer value of relation ( TTOb 
computed with the median value of the wavelength in each in- 
terval. The remarkable agreement demonstrates that the con- 
cept of effective independent coupling channels is physically 
relevant to describe continuously distributed losses. Note that 
the number of channels can also be obtained through the first 
two moments of the width distribution ©: 



M = 



(i 2 )l(iY 



(ii) 



worst case (M = 3) the discrepancy is only of 17%, consistent 
with the relative variation of M given above. 

Thus, having validated the description of losses through the 
introduction of effective scattering channels, whose number 
M is solely fixed through the length of the absorbing part of 
the boundary, the numerical distribution of the complexness 
parameter X is finally compared to the theoretical prediction 
(Fig. |3j for the same absorbing lengths, conductivities, and 
frequency ranges as those given in Fig. [2] The excellent agree- 
ment, even in the tail of the distribution, confirms that the 
prediction Q, obtained within the perturbation theory, con- 
tains the essential features to account for the complexness of 
the resonance states due to spatially continuously distributed 
losses. 

In conclusion, we investigated the statistics of complex res- 
onance states in an open 2D chaotic microwave cavity by solv- 
ing the Maxwell equations in a cavity whose boundary is sub- 
ject to Ohmic losses. We successfully compared the statistics 
of its resonance states and associated widths with the predic- 
tions of an effective non-Hermitian Hamiltonian model. In the 
limit of weak coupling, we have shown that spatially continu- 
ously distributed losses could be mapped to a discrete model 
involving a finite number M of coupling channels, which con- 
stitutes a variable parameter in the cavity. To our knowledge, 
these results are the first unambiguous confirmation of the ad- 
equacy of an effective Random Matrix description to account 
for the spectral and spatial statistics of an open chaotic wave 
system with continuously distributed losses. 

We wish to thank Laurent Labonte for his helpful support 
with Comsol™. 



We checked that except for the smallest value M = 3 the 
estimates ( TTOb and (fTTb agree within a few percents. In the 
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(a) resonances 300 to 400. 



(a) resonances 300 to 400. 




(b) resonances 300 to 400. 




(c) resonances 700 to 800. 
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(d) resonances 700 to 800. 
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(c) resonances 700 to 800. 
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FIG. 2. (Color on line). Distributions of the widths 2 = 7/ {7) 
for different absorbing lengths, conductivities and frequency ranges: 
(a) l abs = 7r/18, cr c = 80Sm _1 and from the 300th to the 400th 
resonances; (b) labs = it/ 6, a c = 400SITT 1 and from the 300th 
to the 400th resonances; (c) l a \, s = tt/6, a c — 400Sm~ 1 and from 
the 700th to the 800th resonances; (d) l abs = n/2, a c = lOOOSm" 1 
and from the 700th to the 800th resonances. The histograms show the 
numerical distribution. The solid lines correspond to the Xm l aw ©• 
The number M of channels used for the comparison corresponds to 
the nearest integer value of relation dlOt computed with the median 
value of the wavelength in each frequency interval. 



FIG. 3. (Color on line). Distributions of the complexness parameter 
Vff E (X) for the same absorbing length, conductivities, frequency 
ranges as given in Fig. [2] The histograms show the numerical dis- 
tributions. The solid lines correspond to the theoretical prediction 
Jilt where the number M of channels corresponds to the nearest in- 
teger value of relation < | X 0| > computed with the median value of the 
wavelength in each frequency interval. 



